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TECHNICAL NOTE NO . 1530 


PLASTIC BUCKLING- OF A RECTANGULAR PIATE 
UNDER EDGE THRUSTS 
By G. H. Handel man and V. Prager 

SUMMARY 


The fundamental equations for the plastic buckling of a rectangular 
plate under edge thrusts are developed on the basis of a new set of 
stress -strain relations for the behavior of a metal in the plastic range. 
These relations are derived for buckling from a state of uniform 
compression. The fundamental equation for the buckling of a simply 
compressed plate together -with typical boundary conditions is then 
developed and the results are applied to calculating the buckling - loads’ 
of a thin strip, a simply supported plate, and a cruciform section. 
Comparisons with the theories of Timoshenko and Ilyushin are made. 
Finally, an energy method is given which can be used for finding 
approximate values of the critical load. 


INTRODUCTION 


This paper is concerned with the plastic buckling of a rectangular 
plate which, previous to buckling, Is under a uniform compressive 
stress ct 0 in the direction of one of its edges . In the case of elastic 

buckling, in which a 0 remains below the elastic limit of the plate 

material, it Is well known that the buckling stress depends on the 
dimensions of the plate and on the manner in which It is supported 
(cf . reference 1, ch. j ) . In the case of the plastic buckling of beams, 
on the other hand, Engesser (reference 2) and Yon Karman (reference 3) 
developed a satisfactory theory based on the fact that for a fiber 
which is compressed beyond the elastic limit the tangent modulus 
(i.e., the ratio of the variation of strain to the corresponding 
variation of stress) assumes different values depending on whether the 
variation of stress constitutes an increase or a relief of the existing 
compressive stress . 

Generalization of this theory to the plastic buckling of plates 
has repeatedly been attempted. These attempts can be divided into two 
groups which may be labeled formal and analytical generalizations. The 
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formal generalizations start from the remark that the formulas of the 
Engesser-Yon TTa-rWri theory of the plastic buckling of beams differ from 
the -we 11 -known formulas for the elastic buckling of beams only by the 
fact that the so-called "reduced modulus” replaces Young's modulus. A 
formal generalization of the Engesser-Yon Karman theory to the plastic 
buc kling of plates is therefore obtained by introducing the reduced 
modulus into the formulas for the elastic buckling of plates in such a 
manner that the results of the Engesser-Yon Karman theory are obtained 
in the case of a narrow rectangular strip which is free on its long 
edges and simply supported on the short edges where it carries a 
compressive load. Of course, this formal generalization is more or less 
arbitrary and leads by no means to a unique result . Formulas of this 
type have been suggested by Bleich (reference k, p. 216 ff .) and 
Timoshenko (reference 1, p. 3®^) • 

In contrast with these formal generalizations of the Engesser- 
Yon Karman theory, the analytical generalizations do not merely introduce 
the reduced modulus of the theory of beams Into the formulas for the 
elastic buckling of plates. Instead, the analytical generalizations go 
back to the considerations by which the reduced modulus Is derived and 
try to apply these to the case of a buckled plate. Generalizations of 
this kind have been previously presented by Eaufmann (reference 5) and 
Ilyushin (reference 6) - As is shown in the present report, however, 
these authors use stress-strain relations which do not fulfill certain 
postulates of the theory of plasticity; the correctness of their results 
must therefore be questioned. 

The present paper aims at developing a theory of the plastic 
buckling of plates which takes full account of the modern theory of 
plasticity. The stress-strain relations in the plastic range are 
discussed at considerable length in the first section of the ANALYSIS, 
and it is shown that, for an adequate treatment of buckling phenomena, 
a theory of plastic flow Is indicated rather than a theory of plastic 
deformation of the type used by Eaufmann and Ilyushin. The precise 
definitions of these terms and the basic considerations suggesting the 
use of a theory of plastic deformation for problems such as buckling are 
fully discussed in the ANALYSIS.' A particular theory of plastic flow 
suitable for the treatment of the problems under consideration is 
developed in the first section and Its relations with other theories of 
plasticity are pointed out. It Is shown that in the particular case of 
a plate buc kling out of a state of simple compression there is very 
little freedom in the choice of the Btress -strain relation if It Is to 
fulfill certain simple postulates. This means that all the empirical 
Information which is necessary for the theoretical treatment of the 
plastic buckling of a rectangular plate under edge thrusts can be 
obtained by a simple compression test. 

The second section presents the development of the fundamental 
equation of the plastic buckling of a simply compressed plate, and the 
appropriate equations describing typical boundary conditions are given 
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in the third section. The remaining parts contain several e x a mp les, 
■which are carried out in detail, as well as an equivalent energy 
principle which proves to he very useful for approximate computations • 
Finally, the appendixes contain detailed discussions of several 
technical points raised in earlier parts of the paper. 

The authors are indebted to Professor H. S. Tsien for drawing their 
attention to the work of Kaufmann (reference 5) * 

This investigation was conducted at Brown University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SIMBOIS 


a 

A 

A 


length of plate 

constant in compressive stress -strain law- 
expression in variational principle 


a', a", a' ", b*, b", b ,M , c‘ 


coefficients in plastic 'stress -strain 
law 


b 

B 


width of plate 

expression in variational principle 


c = a - l )/[(5 - hv)\ - (1 - 2V)5 


c l, 62, eg, c^ arbitrary constants appearing in equation for w 
D flexural 


rigidity of plate (#*=/(?• - v2 )) 


D 


11' D 12 j B 2Q 


coefficients in plate equation for plastic flow 


D ll' = d 12' = ^ 22 . 1 ^ > D 22 I = ®22/ D 

P-jj*, D 1P *, D22* coefficients in plate equation for Ilyushin's theory 

of plastic deformation 

D 11 = D ll* Z 15 ' D 12 = ■ D 12* / ■ D 22 = 

E tangent modulus in compression 
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E 


o 


E* 


Es 
t( y) 

h 


I 


Young 's modulus 

Yon Karman's reduced modulus 

secant modulus obtained from compressive stress -strain diagram 
section of buckled middle surface for x = Constant 
thickness of plate 
moment of inertia of cross section 


«o a 


k =12 


k - ^k 
k - a 2 


! (l - V*) 


* 2 h% 0 


& = (2 - Vji-L + (2V - l)Kg 

a a^/ax 2 , Kg = a 2 v/a y 2 , ±22 


2a 2 w/a*a y 


m 


M-v 




M 


^7 


N. 


N-« 


N t 


number of half waves in buckled configuration 
rate of change of bending moment about y-axis 
rate of change of bending moment about x~axis 
rate of change of twisting moment 
integer 

reduced compressive stress resultant 

rate of change of stress resultant in x-direction 

rate of change of stress resultant in y-direction 


p 2 = (m) 2 Biz 

D 22 


P 


total compressive force (c 0 bh) 
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Q = [Dig - (1 - Vttj/Dgj, 

2 = fo^ _ fm^ 2 D 11 :d 22 ~ p 12 2 
r = ^p(q. + p) 

R side ratio (b/a) 

e = \/p(q. “ p) 

t tine 

u = 8V 2 + 12v - 23 

v deflection rate 

i, j, z rectangular Cartesian coordinates; x,y-plane coincides with 

middle surface of unbuckled plate 

(2 - v)e + (2V - l) €p 

z = i £ 

° ± 

a = 1 ' o{5 - 4 V 1 ) 

a. ratio of Yon Harman 'b modulus to Young's modulus (in section 

"Buckling of a Bimply supported plate" only) 

p constant in compressive stress -strain law 

± >0 
£ <0 
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uniaxial strain 


de x> de y> &€ z’ 

d ^yz* d ^: 


infinitesimal strain increments present in buc kling 


zx 


de x’> V 

' d€ z 

' dy ' 

d 6 s"' d* y 

de z 

", d 7 " 

’ 'xy 

« x (3:,y) 

normal strain 

« 2 ^ x >y) 

normal strain 

?(x,y) 

shear 

strain 2 

i o - Sz o/ h 



to + - “ * ' 

\Ja? - 

1 

£ ss -a - 


• 1 


tj = sb 

K 

X, 

V 

I = rb 
a 


ratio of Von Karman 's modulus to Young's modulus (E*/E 0 ) 
ratio of Young's modulus to tangent modulus (E^E) 
Poisson 's ratio 


uniaxial stress 
intensity of stress 


(fJT 


a x?y 


+ 


CT cr critical compressive stress 

<y Q = -d x original compressive stress in plate 

a x , Oy, a z normal stress components 
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Tjy, T yz’ T zi shear stress components 

CLCf-y , d0y ^ J 

Infinitesimal stress increments present in buckling 



cb function of the intensity of stress 

cd‘ = dco/dcr^ 
fl = oj/(l - cd) 



SI* = d^dn 1 

superscript denoting values on unloading side of neutral 
surface 

+ superscript denoting values on loading side of neutral surface 

ANALYSIS 

Stress-Strain Relations for Buckling from, a 
State of Uniform Compression 


The mechanism of buckling beyond the elastic limit is relatively 
complicated because the material, -which was originally in a state of 
simple compression, is loaded in some regions and unloaded in others 
during the buckling process. Consequently, the stress -strain relations 
must be considered in some detail with special reference to the problem 
of loading beyond the elastic limit followed by unloading. 

The material must exhibit strain-hardening if the determination of 
the buckling stress is to constitute a problem. Indeed, for a perfectly 
plastic material which yields under constant stress. Yon Kerman 's reduced 
modulus vanishes once the initial compressive stress has reached the 
yield limit. This means that the bending stiffness is reduced to zero 
and buckling must be expected quite independent of the dimensions of the 
bar. 
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Stress-strain laws for materials which exhibit strain -hardening 
can he divided into two types which, for convenience, will be called 
"theories of plastic deformation" and "theories of plastic flow*" 
According to the first group, there exists a one-to-one correspondence 
between stress and strain in the plastic range, as well as the elastic, 
provided that the material is being loaded. The stress-strain law of 
the well-known Henciy-Nadai theory (reference J, ch. li)-, and reference 8) 
and the law used by Ilyushin (reference 6) in his discussion of plastic 
buckling are typical theories of plastic deformation. On the other 
hand, the theories of plastio flow are based on the assumption that, 
for a given state of stress, there exists a one-to-one correspondence 
between the rates of change of stress and strain in such a manner that 
the resulting relation between stress and strain cannot be integrated 
so as to yield a relation between stress and strain alone. Typical 
examples of theories of plastic flow are the stresB-strain relations 
developed by Prager (reference 9 ) and Hands Iman, Lin, and Prager 
(reference 10 ) . A particularly important difference between these two 
basic theories of plasticity lies in the faot that according to the 
theory of plastic deformation the strain which corresponds to a certain 
state of stress the theory of plastic deformation is entirely independent 
of the manner in which this state of stress has been reached, whereas, 
according to the theory of plastic flow, the Btrain depends on the manner 
in which the state of Stress is built up. 

The stress-strain relations to be uBed in the analysis of the 
plastic buckling of a rectangular plate under edge thrusts form a 
special case of those developed by Hands Iman, Lin, and Prager in 
reference 10. In this particular case, however, it is possible to 
develop the stress-strain relation in a quite elementary manner, an<£ 
the inherent difficulties of the theories of plastic deformation can 
be seen from a slightly different point of view. It appears worth 
while, then, to examine these relations in some detail with special 
reference to the problem which forms the subject of the present report. 

In the following, the stresses and strains in the buckled plate 
will be referred to a fixed system of rectangular Cartesian coordi- 
nates x, y, and z. The x,y-plane of this coordinate system coincides 
with the middle surface of the unbuckled plate, and the axes of x 
and y coincide with two of its edges, the other edges falling on the 
lines x = a and y = b. Prior to buckling, the plate is under a 
uniform compresBive stress a 0 in the direction of the x-axis (fig. 1) . 

The investigation of the stability of the state of stress 


a 


x 


-a 


o 


a 7 " 


a z = o 


sy 


= T 


yz 


= T 


ZX 


= 0 



« 

* * 


( 1 ) 
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requires the knowledge of the relations between the Infinites imal 
Increments of stress dcr s , dcr y , a a z , dT.^, dT yz , and dT zz and 

the corresponding increments of strain de„, de_, de„, d 7 v _ r , d7_„, 
and d 7 zx - Within the framework of plate theory, cr z = = t^. z = 0 , 

even in the buckled state, and hence dcr z = dT Y7 = dT^ z = 0. Accord- 
ingly, d7 xz = d 7y Z = 0. Within the elastic range the remaining 
increments of stress and strain are related to each other by means of 


E 0 d « x = dcJ x " V ^y 

E o d€ y = " V + ^y 
E q de z = -V db x - V dn y 

Eq d 7 = 2(l + V) dT-jjy 


( 2 ) 


idler e E Q denotes Young's modulus and v , Poisson's ratio. Before an 

analysis of the plastic buckling of the plate can be attempted, the 
relations replacing equations (2) in the plastic range must be known. 

In order to establish these relations, it will be convenient to think 
of the strain increments as consisting of reversible (elastic) and 
permanent (plastic) components : 


d *x = d6 x' + d€ x" 


d6 y “ d 6 y ? + d 6 y' 


d« z = de z ' + dV 


( 3 ) 


d V - d V + d V 


Primes and double primes denote elastic and plastic conqponents, 
respectively. The elastic increments of strain are related to the 
increments of stress by means of equations (2), in -which the left-hand 
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E 0 de x ' = dff x - V a 0 y 
E 0 d6y ' = -V do x + da y 
E 0 de z ’ = -v dcr x - v dc y 
E o «v = 2(1 + v > 


W 


These relations may he regarded as the definitions of the elastic 
increments of strain- The purpose of the following discussion is to 
establish similar relations for the plastic increments of strain. 

The elastic increments of strain, equations (!<•) , depend only on 
the increments of stress and are independent of the existing stress cr 0 

Moreover, a reversal of the signs of all increments of stress leads to 
a mere reversal of the signs of all elastic increments of strain. The 
plastic increments of strain, however, do not have these properties; 
since they must vanish as long as a Q remains below the elastic 

limit, they cannot be independent of the existing stress a Q . 

Moreover, if for a given value of a Q certain stress increments 

produce plastic increments of strain, stress increments of the same 
magnitudes but opposite signs do not produce any plastic deformation. 

In other terms, beyond the limit of elasticity an infinitesimal change 
of stress may be classified as loading the material or nob according 
to whether it is accompanied by permanent deformation. Infinitesimal 
changes of stress which do not load the material may be classified 
in turn as unloading or neutral. Unloading brings the material into 
a state of stress such that all sufficiently small further changes of 
stress are accompanied by elastic deformations only. These basic 
differences in loading and unloading appear somewhat more clearly if 
the simple example of a uniaxial state of stress and strain (say a 
tensile test) is considered. Let a denote the stress and e the 
strain in figure 2 and suppose the material is loaded to the point P. 
The stress-strain diagram for unloading is a straight line PA with 
the same slope as the loading curve at the origin 0. The permanent 
strain corresponding to loading up to the point P is measured 
by 0A . Suppose now that, after the point P has been reached, the 
test specimen is further loaded to the point P^j by this, the 

permanent strain is increased by the amount AA-j_. In other words, the 
change from P to P^ constitutes loading in the sense just defined. 

On the other hand, if the new state of stress and strain is given by 
the point Pg, that is, if the stress has been reduced below that 
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at P, the permanent strain is left unchanged- Furthermore, any 
small change of stress from the point Pg (strictly, all changes 

■within the ranges PrP and PgA) produces an additional deformation 
•which is purely elastic. The material has thus been unloaded. 

For uniaxial stress any change of stress constitutes either 
loading or unloading. A third possibility, designated as 'Watral" 
change of stress, exists in the case of combined stress. A neutral 
change of stress, -while not accompanied by a permanent deformation, 
brings the material into a state such that there exist certain further 
changes of stress •which are arbitrarily small and yet produce a 
permanent deformation. This third condition is illustrated in the 
analysis of the buckling of a plate. It. is precisely the possibility 
of the occurrence of neutral changes of stress -which distinguishes 
the present problem from that treated by Engesser and Ton Ka'rWun, 
for in their case the stress is uniaxial in the buckled state as 
veil as in the unbuckled. Accordingly, a change of stress can be 
only an increase of the existing compressive stress (loading) or a 
decrease (unloading) . The situation Is more complicated in the case 
of a plate. 

Since there Is no permanent deformation accompanying neutral 
changes of stress or unloading, d6 x " = de y " = de z " = d7™. " =0; and 

the relations of equations (4) define the total change of strain. For 
loading, hovever, equations (4) must be supplemented by equations of 
the form 


E„ 

d€ x " = a ' d0 I + b * 

do 

o 


y 

®o 

d« y " = a” dff z + b" 

ddy 

E_ 

ds " = a " ' dcr + b 

”'d0. 

o 

z X 


E o 

$ 

ii 

IO 

o 

9 

3 



(5) 


•where the coefficients a', b', a", b", a'", b"', and c' depend 

on the existing stress a Q . 

As is customary in the theory of plasticity, the plastic defor- 
mations vill be supposed to represent a mere change in shape but no 
change in volume. Accordingly, 



12 


NACA TN No . 1530 


d 6x " + d6 y " + de z " = 0 ( 6 ) 

This relation must hold independently of the values of dc z and defy • 
Thus, 

a' + a" + a'" = 0 (7) 

and 

b' + h 11 + b'" = 0 (8) 


The elastic formulas, equations (4), exhibit a certain symmetry of 
the coefficients appearing on the right-hand side. For instance, the 
coefficients of dcr^. in the second and third equations are equal, as 

are the coefficients of day in the first and d 0 z in the second 

equations • Which of these symmetries, if any, will be maintained in 
equations (5) ? The existing state of stress singles out the x-axis, 
but it does not matter which of the other two axes is labeled y 
and which z . Accordingly, 


(9) 


In view of equations ( 7 ) and (9)> 


a 


ff 


ft r 

a 



(10) 


These coefficients can easily be expressed in terms of the so'- called 
"tangent modulus" corresponding to the compressive stress cr 0 . 

Application of equations (4) and (5) to simple compression in the 
x-dlrection yields (with defy = 0) 


E o + a*/) 

= dcr z + a 1 ckr z 


( 11 ) 


or 


to x E o 

d 6 x = 1 + a’ ■ 


E 


( 12 ) 
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■where E denotes the tangent modulus. With 


X = E 0/ /l 


equation (ll) gives 


a' = X - 1 


( 13 ) 

(HO 


Hence, according to equation (10), 


a " = a”’ = 


(15) 


Next, the criterion for neutral changes of stress must be considered. 
Any given Infinites ima l change of strain can be decomposed in the 
following manner: 


de x = -i(de x + de y + de z ) + ^(2de x - de y - de z ) 
de 7 = K d6 x + d€ y + d€ z) + K" de x + 2de y " 


d6 z - i( d6 x + d6 y + de z) + K' 


de~ - de^ + 


2dej) 


d7 sy “ 0 + d7 xy 


d7 yz — 0 + d7 yz 


d7 zx = 0 + dr 


zx 


(16) 


The change of strain defined by the first members of the right-hand 
sides of these equations is a uniform expansion (or contraction) in _ 
all directions . Such a uniform expansion changes the volume but not 
the shape of the element to which it is applied. The change of strain 
defined by the second members, of the right sides of equations (l 6), on 
the other hand, affects the shape of the element but preserves its 
volume . The work done by the existing stress or stresses on the change 
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of strain, equations ( 16 ) , consists of the work done on the change of 
volume represented by the first members of the right sides of 
equations (l6) and the work done on the change of shape represented 
by the second members . Since all changes of volume are supposed to 
be of an elastic nature, it seems natural to speak of loading or 
unloading according to whether the work dW which the existing 
stresses do on the change of shape alone is positive or negative. 
Vanishing of this work must then be interpreted as indicating a 
neutral change . 

In the case under discussion, the only existing stress 
is a x = ~ a o> the criterion for loading or unloading is furnished 

by the sign of 


£W 


- _ (2de x - d - de z ^ 


(17) 


■while neutral changes are characterized by 


24 6 X - d€y - de z * 0 


( 18 ) 


Now, for unloading, the entire change of strain is of an elastic 
nature and equations (2) apply. Equation (17) is therefore equivalent 
to 

E 0 dW = - v dxr y ) - (-v d 0 x + d by) - (-V da x - V d 0 y )J (19) 

Since cr 0 >0, this expression will be negative, whenever 

2 d 0 x - dcr y >0 (20) 


This inequality, equation (20), is thus seen to constitute the criterion 
for ■unloading. Similarly, the criterion for neutral changes of stress 
is found to be 


2dcfx “ 


dhy = 0 


( 21 ) 
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Changes of stress ■which satisfy neither equation (20) nor (21), that' 
is, changes of stress for which 


2dn x - doy < 0 


( 22 ) 


must therefore constitute loading. Another definition for the criterion 
for the three types of change of stress, which is found "by combining 
equations (20), (21), and (22), is that the change of stress is classified 
by the sign of the increment in the second invariant of the stress 
deviator, which measures the intensity of stress • A detailed account 
of this alternative formulation is found in reference 10 . 

By a suitable choice of the values of day and dcy, the 
expression 2da x - doy can be made . to fulfill the following 
inequalities : 


0 > 2dcr x - day > - e 


-(S3) 


where g is an arbitrarily prescribed small positive number. All 
changes of stress satisfying equation ( 23 ) constitute loading and are 
therefore accompanied by plastic deformations in accordance with 
equations ( 5 ) • For e — > 0, however, these changes tend toward 
neutral changes of stress for which there are no plastic deformations. 
Furthermore, there Eire no plastic increments of strain when 
2da x - day > 0 . It is to be expected that the total strain increments 

will be continuous in the region which marks the transition from 
unloading through, the neutral state to loading. Although such a 
statement does not follow specifically from the equations of equilibrium 
or compatibility, continuity should be expected in the strain increments . 
With this assumption, the plastic increments of strain, equations ( 5 ), 
should vanish whenever the increments of stress satisfy equation (2l) . 
This furnishes the conditions 


a' + 2b' = 0 
a" + 2b" = 0 
a m + 2b'" = 0 
o' = 0 


( 21 ) 
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Together with equations (lb) and ( 15 ) , these equations determine all 
coefficients appearing in equations ( 5 ), which therefore take the 
form 


E 

E 

E 

E 


o 


o 


o 


o 


d6 x " = (X, - 1) dcr x - - g - dc y 

ae z" - + Hr 1 

d V'° 


(25) 


It is interesting to note that here again the coefficients of do x in 

the second and third equations are equal, as are the coefficients 
of dcr y in the first and of &a x in the second equation. Whereas 

in the elastic case this type of symmetry in the stress-strain relations 
is a consequence of the isotropy of the material, this is no longer so 
in the case of equations (25)* Indeed, the equality of a" and a'" 
(see equation ( 15 ) ) follows from the assumption that the plastic 
deformations do not involve a change in volume . The equality of b 1 
and a", on the other hand, might be described as almost accidental, 
the value of the ratio a "/a 1 being fixed by the assumption Just 
mentioned, while the value of the ratio b '/a ' is fixed by the form 
of the condition for neutral change of stress. 

Combination of equations (4) and (25) finally yields the stress - 
strain relations which will be used throughout this paper: 


1 , 


d e x - * ^x " ( 


d v ’-( V + H 

X ) <^x + HP to ; 

d «z- -( V + H' 

-) aox - ( v - H 

«L 7^ = 2(1 + V) 

d V 


(26) 


'y 
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It seems -worth, -while to stress once again the assumptions on the 
basis of -which these stress-strain relations are derived. These are 

(1) Plastic deformations do not involve a change of volume 

(2) The criterion for loading or unloading is furnished by the 
sigp. of the -work dW -which the existing stresses do on the change 
of shape produced by the increments of stress 

The first assumption is conmonly made in the theory of plasticity 
(cf . reference J, p. 10} and is confirmed by the exper im ents of 
Bridgman (reference 11, p. 1 66 ). The second assumption is a slight 
generalization of a similar assumption which Prager (reference 10) 
introduced in the case of incompressible plastic materials; more 
recently, it has been used by Ilyushin (reference 6) . 

It is interesting to note how far the stress-strain relations, 
equations ( 26 ), differ from those used in previous work on the plastic 
buckling of plates. In the present notation, Kaufmann's stress-strain 
relations (reference 5 ) are 


E de = \ dcr - V fay 
o x x y 

E q dgy = -XV dcr x + dcr y 
E q dr^ = (1 + X)(l + V) dr^ 


(27) 


(The expression for d€ z is not given because this strain component 

is nob necessary for the determination of the bending and twisting 
moments in the buckled plate . ) . It is seen that here the coefficients 
of dcr y in the first equation and that of dn x in the second are 

unequal. In an earlier paper on the plastic buckling of cylindrical 
shells (reference 12, footnote 1, p. 422) in which similar stress- 
strain relations were used, Kaufmans comments on this lack of symmetry, 
recommending that the stress -strain relations, equations (25), be 
checked by experiment. Since this type of symmetry in the present 
stress-strain relations, equations ( 26 ), has been characterized as 
almost accidental, the lack of symmetry in Kaufmann's relations hardly 
constitutes a sufficient reason for discarding the stress-strain 
relations, equations (27) * It is not difficult, however, to show 
that these relations correspond to an unacceptable condition for 
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neutral changes of stress. Indeed, subtraction of the elastic 
increments of strain, equations (k) , from the total increments of 
strain, equations (27), yields the following plastic increments of 
strain: 


E q de z " . (X - 1) 6a x 
E q de y " = v(X - 1) dc x 

E o d7 *y " “ ft “ + v ) dT xy 


(28) 


These plastic increments of strain vanish if 


d0 x = 0 


dT 


37 


= 0 


n 


(29) 


According to Kaufmann’s stress-strain relation, neutral changes of 
stress are characterized by the two conditions given as equations (29) • 
If the most general change of stress considered here is represented 
by a point with the coordinates dc^, dcr y , and dr ^ in a three- 

dimensional space, the condition of equation (21) represents a plane 
through the origin which separates the "region of loading" from the 
"region of Tinloading." Equations (29), however, define a straight 
line which does not mark off two such regions. 

Ilyushin (reference 6) considers an incompressible material and 
assumes the stress-strain relations for loading to have the form 


s = 


2cr„ - cr 


° x 2(1 - cp 


f 


E„e„ = 


2a y - a x 


o y " 2(1 - a) 


( 30 ) 
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■where cu is a function of the intensity of stress defined "by 



For loading, the increments of stress and strain are then connected "by 

1 f 2n x — ff y 

E^ dg_ = ■ryt r{2dc_ - dto - — den) 

o tLfc x 2(1 - ai) x ^y 1 - oj J 

E ° 4 v - ‘ to ) 

E o - r^v - r?z *°) 

For buckling from, a state of uniform compression cr x = -a Q , in 
particular, 

to = d0 ± = -^(2dcr x - dc y ) (33) 

•where 

cd ' = to/drr^ (3*0 

Equations (32) then reduce to 
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For unloading, the relations, equations (2), are supposed to hold with 
V = l/2 on account of the assumed incompressibility of the plate 
material: 


E o de x = i( 2da x “ ^y) 

E ° d< V = 2( 2d0 7 ‘ 

®o d ^xy = 3dTjy 



(36) 


As to the criterion for loading and unloading, this is again supposed 
to he given by the sign of the expression 

dW = a x d€ x + a y ds y + r ^ d 7 ^ (37) 

In particular, it is given by the sign of 

-cr 0 de x (38) 


in the case of buckling from a state of uniform compression cr x = -cf Q . 

In view of the first of equations (35)> this means that neutral changes 
of stress are again characterized by equation (21) . It is easily seen, 
however, that for 2da x - defy = 0, equations (35) &&d ( 36 ) do not give 

the same increments of strain. Ilyushin's stress-strain relations are 
thus seen to exhibit an objectionable discontinuity along the 
surface 2dh x - ddy = 0 which separates the region of loading from 
the region of unloading. A more detailed analysis of the effect of 
this discontinuity in the case of a buckling plate is found in 
appendix A. 
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Fundamental Equation of Plastic Buckling of a 
Simply Compressed Plate 

The technique used in the derivation of the fundamental equation 
of the plastic buckling of a simply compressed plate is quite similar 
to that needed for the same problem in the elastic range- (See, for 
example, the more general problem of combined bending and compression 
of elastic plates in reference 1 , p- 302 - ) There is one essential 
difference, however, in that the stress-strain relations given in 
equations (26) must be used in the regions of loading rather than 
generalized Hooke 's law- Consequently, the middle plane of the 
unbent plate will no longer play the role of the neutral surface in 
the buckled position. Once the position of the neutral surface has 
been found and the bending and twisting moments determined as functions 
of the second derivatives of the deflection of the plate, the 
equilibrium conditions and the final differential equation can be 
derived in exactly the same fashion as that used by Timoshenko in 
reference 1- 


It will be found more convenient, in the following discussion, 
to use "reduced stresses" rather than actual stresses, that is, stresses 
reduced by dividing the actual stress by Young's modulus E Q . No new 

notation will be employed to denote these reduced stresses; therefore, 
care must be taken in interpreting the results obtained here in terms 
of the known facts for elastic buckling. An attempt will be made at 
such points to keep the notation clear. In addition, the use of 
differentials of stress and strain may lead to Borne confusion in 
deriving the equations of equilibrium for an element. Since the 
stress-strain relations gLven in equations (k) and (26) are linear 
in these differentials, both sides of the equations may be divided 
by dt > 0, where t may be regarded as the time. It should be 
noted that t appears homogeneously; that Is, the time scale may be 
arbitrarily distorted without changing the equations . If differenti- 
ation with respect to t Is denoted by a dot, equations (k) and (26) 
can be rewritten as reduced stress-strain relations, for W = 0, 
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and for W = 0, 


where W = dW/dt. 


6 X = di - Va y 


V = ' V<5r x + *y 


€z = - V o x _ va y 


7^ = 2(1 + v)^ 






The stress rates cr x and cr y can he found in terms of the 

corresponding strain rates for loading by solving the first two of 
equations (39) • Thus, 


a _ 

(5 - 4v)X - (l - 2V) 


■ (1 - 2V )g [ (x ♦ 3, 'x * - 1 * £ v) q 




(5 - 4v)X - (1 - 2V) 


-( 2 a - 1 + 2v)e x + bxQ 

- (1 - 2V) 2 L y J 


(to) 


The criterion for loading 


2ct x - cr y <0 


(to) 


can then he written as 


(5 - 4v)X - (1 - 2V) 


? f(2 - v )* x + (2v - Del < 0 

- (1 - 2v) 2 L x yJ 


( 43 ) 


N’ow, Poisson's ratio v satisfies the inequality -1 V = ^ 

(of. reference 13, P- 104) ; in addition, X = 1- Consequently, the 
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expression appearing outside the brackets is always positive and the 
inequality, equation (43), can be replaced by 


(2 - V)e x + (2V - l)e y < 0 


(44) 


The strain rates appearing in equations (39) and (40) must now 
be evaluated- The strain rates in the middle surface -will be denoted 
by 4^ = 4^(x,y), the normal strain rate in the x-direction; 

= e^(x,y), the normal strain rate in the y-directionj and 

y = 7 (x,y), the rate of shear strain- Points on the normal to the 
undeformed middle surface are assumed to remain on the normal of the 
bent middle surface. This implies that the strain rates e 4 , 

and at any point of the plate can be written in the following 

form: 


' ^12 


1 


(45) 


The quantities E-^, Kg, and K-|p appearing in equations (45) are 

defined in terms of the rate of deflection w = v{x,y) of the middle 
surface in the following way: 


-.2- 

K , = £-2 
1 3x2 


Eg = 


Sy 2 


Eig = 2 


-v2- 

o w 


( 46 ) 
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* > 

Geometrically, K]_ and represent the rates- of curvature of 

the middle surface in the x- and y-directions, respectively, 

■whereas Kj2 represents the rate of relative twist. The criterion 

for loading, equation (bb) , can now he rewritten in terms of the 
strains of the middle surface and the quantities K^, K^, and K-|p . 

It is seen that loading takes place provided 

(2 - V)e x + (2V - l)e 2 < z |(2 - V)K X + (2V - 1 )kJ (47) 


With the word "sign" to denote the sign of the quantity within the 
parenthesis and vertical bars to denote the absolute value of the 
enclosed expression, the last inequality may be transformed into 



V)±! 


+ (2V - 1)kJ 


> 


(2 - y)6 L + (2V - l)6 g 
| (2 - V)K X + (2V - l)Kg| 


(48) 


This Inequality can be simplified further by Introducing two new 
quantities E and z Q defined by 


K = (2 - V)± 1 + (2V - 1)Kq 

(2 - V)e 1 + (2v - l)e 2 

z = =~ 

° K 



(49) 


The inequality, equation (48), becomes then 


sign (K) >Ez q J\e \ = z 0 sign (E) 


z 


(50) 
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i 


that is. 


z > z^ 
o 


for positive E 


z < z 

o 


for negative 


E 



(5D 


The surface z = z Q separates the regions of loading and unloading in 

the plate; a given part is in a state of loading or not according to 
■which condition of equations ( 51 ) Is satisfied. 

The criterion Just developed must now be applied to the problem 
of buckling. As mentioned previously, the stress distribution of the 
buckled plate differs from, the original state of pure compression by 

certain additional stresses cf x dt, dy dt, and T^y dt . These new 

stresses are such that their total stress resultants must vanish and 
the moments produced will be in equilibrium with the moment generated 
by the original compressive force in the buckled plate. The vanishing 
of the stress resultants will lead to a formula for z D in terms of 

the constants of the material and the value of cr 0 . Once this equation 

has been developed, a rather straightforward computation will lead to 
the desired equation of equilibrium. 

> • 

The rates N x and Ny of the stress resultants are defined as 


N x = 


*7 = 


>h/2 . 
f-h/2 

>h/2 
1 -h/2 


q x dz 


dy dZ 


(52) 




As indicated in equations (51)? two cases must be considered according 

to whether E > 0 or E < 0 . 3Tor E > 0, direct computation shows 
that 
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N x - vN y = «ih + cK^| - z 0 ) 2 
% - VN X = e 2 h - icK^| - z Q ) S 





(53) 


where the quantity c is a function of X and V given hy 


c 


X - 1 

(5 - 4v)\ - (1 - 2v) 2 


( 54 ) 


Appendix B contains the details of this calculation and others used 

in this section. It has been pointed out previously that N_ and N, 

A « 

must vanish. According to equations (53) , this yields 



(55) 


Thus the strain rates 4 ^ and in the middle surface are 

related by the equation 

(56) 

From the definition of z Q , equations (49) , and this result, it is 
seen that 



= (2 - V)-^ + (2V - 




4v) 


(57) 
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This result may then, he substituted hack into equation (56) to yield 




2g gh 

ct 


c (5 - i 4 -v) 


(58) 


Relation (58) can he solved for z Q to yield the equation of the 

neutral surface ■which separates the regions of loading and unloading. 
Since h, c, and V depend only on the geometry of the plate and 
the applied compressive stress, z Q will depend only on these 

quantities. It is more convenient to introduce a new quantity 
defined hy 



.( 59 ) 


Then the quadratic equation for z Q , equation (58), becomes a 
quadratic equation in £ Q , namely. 


£ 



+ 


c (5 - ^v) 


0 


(60) 


There are two solutions to this equation, in general, hut the only 
one which is physically realizable is 


£0 = £ o + =• a + \/® 2 " 1 


(61) 


where 


a = 1 


2 

0(5 - 4 -v) 


(62) 


Equation ( 6 l) gives the desired formula for the neutral surface. 



28 


NACA TN No . 1530 


When E < 0, the procedure is exactly the same as that outlined. 
Again the details are found in appendix B. The formula for the neutral 
surface is given in this case by 




(63) 


9 

Roughly speaking, the sign of E indicates -whether the plate "buckles 
up " or '"buckles dovn . " Consequently, the differences in sign found by 
comparing equations ( 6 l) and ( 63 ) are quite natural. It would also be 
expected that the rates of change of the bending and twisting moments, 
as well as the resulting equilibrium equation, should be independent 
of the sign of E. This will be shown to be true. 

The rates of change of the bending and twisting moments can be 

• # 

computed now that z Q or is known for E > 0 and for E < 0 . 

The rates of change of the bending moments, and My, ■ are 

defined as 


My 



a x z dz 


dyZ dz 


(64) 


where the moments are taken about the y- and x-axes, respectively. 
The rate of change of the twisting moment M_y is given by 





t z dz 
xy 


(65) 


The calculation of the rates, equations (64) and (65), must be carried 
out separately for E> 0 and E < 0 • It can be shown that the only 
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quantity appearing in the final result -which depends on the sign 
of K is the function 8 defined hy 


8 -i 


for K > 0 




8 = 1 
2 

i + |^o - t(fo)~ 

for E <0 


( 66 ) 


According to equation ( 63 ), £ = - £ Q , 

obtained from equations (66) will therefore he the same in either' case. 
Thus the expressions for M^., My, and will he the Bame in both 
cases. The details are found in appendix B in which it is shown that 


the numerical Value of 8 


K ■ - ^ T^) Jjj 1 " 08(2 - V) 9 + i 2& ' c6(2 - V)(2V - 

My - - ~ lg ( 1 ^ 3 , ^2) ‘ c8(2 ■ v) < 2V - l)] + Kg [1 - CS(2V - 

= 12.il- v2) [ (1 2 V) 



( 69 ) 


The equation of equilibrium can be set up in terms of the bending 
and twisting moments and compressive load without reference to the 
stress-strain relations. This has already been done by Timoshenko 
(reference 1, p. 305) for the more general case of combined bending 
and tension or compression. His results may be applied to this special 
case of a simply compressed plate. With the present notation, the 
equation of equilibrium is 
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2 


^xdy + 


S 2 My 

Sy 2 



(70) 


Timoshenko 's relation was originally written in terms of the actual 
bending moments and actual compressive stress resultant N x rather 

than the rates of the reduced quantities. Timoshenko's equation can 
be differentiated with respect to time and divided by S c on both 

sides, so that equation ( 70 ) is the desired equation of equilibrium 
provided N x is defined as 


N x = c 0 h/E 0 (71) 

From equations ( h-6) , 

K ± = 

Kg = d 2 w /dy 2 
K 12 = ^^v/dxdy 


With these relations and equations ( 67 ), (68), and ( 69 ), equation (70) 
may be rewritten as 






( 72 ) 


where the subscripts denote partial differentiation with respect to 
the variable named and 
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Dll = D 

[l - c5(2 - V) 2 ] 

d 12 = 15 

[l - c5(2 - V)(2V 

D22 ~ 

[l - c8(2V - 1)^] 

D = 

.h 3 ! 

0 


12(l - V 2 ) 


(73) 


The quantity D is the well-known flexural rigidity of the, pl^te . 
Equation (72) resembles the equation for the buckling of aif* anisotropic 
plate (reference 1, p. 3^0)' jThere is one important difference, 
however. In the case of ar£ anisotropic plate, the coefficients 

D-|P, and Dpp are constants of the material; for the plastic 
case D D'-jp, and Dpp are functions of cr 0 . In other words, 

the plate is anisotropic but this anisotropy is .caused by and Is a 
function of the compressive stress. Consequently, certain changes 
must be made in the standard procedure for calculating buckling loads 
for anisotropic plates. Several examples Illustrating this technique 
are given in the succeeding sections . Graphs of the quantities 


and 


D n ' = 1 - c8(2 - V) 2 
D ' = 1 - c6(2 -V)(2V - 1) 

D 22 ' = 1 - c8(2 V - l) 2 


as functions of X for v = 0-32 are given in figure 3 and the 
numerical values are listed In table I. 
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Typical Boundary Conditions for the Fundamental Equation 
for a Simply Compressed Plate 

The discussion of the boundary conditions for the buckling 
equation, equation ( 72 ), is facilitated by expressing the moment 
rates M x , M y , and in terms of the second derivatives of 

the deflection rate v and the stiffnesses D-^, ^12> and ^22 

introduced in equations (73) • Thus, it follows from equations ( 67 ), 
( 68) , and ( 69 ) that 


% “ " D lAx “ [?12 “ C 1 - v)D]w yy 

E o*V = -[ D 12 " ^ ' D 22*yy 


( 74 ) 


E Ay - - v )* : 


*y 


The following boundary conditions are typical in the buckling of 
rectangular plates: 

(1) Simply supported edge at x = 0*- The deflection rate w 
and the moment rate must vanish at this edge; that is. 


w = 0 

" D li* r xx ” [?12 ' C 1 - v ) E ]'”yy = 0 


for x = 0 • 

(2) Built-in edge at x = 0.- At -this edge, the deflection 
rate w and the slope rate w x must vanish; that is, 


(75) 


w = 0 
w x = 0 



for x = 0 . 


(76) 
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(3) Free edge at y = 0.- For a free edge, the rates of the tending 
moment My. and of the equivalent shear load ^-£^14 ^ j c5z) + yj^j) 
must vanish (reference 1, p. 300 ) . Consequently, 


[ D 32 - t 1 - V > D ]tac * B 22V * 0 

[ d o£ + (1 ' + “aaV, - ° 


(77) 


for y = 0 . 

(k) Plane of symmetry at y = 0.- If the buckled shape of the 
plate is symmetrical with respect to the plane y = 0, the rates of the 
slope Wy and of the equivalent shear load ^-£3 m 

will vanish. Therefore, -- . 


= 0 


[“02 + ^ - V > D ]* o 7 * D 22 - OT - 0 




( 78 ) 


or 


v y - 0 


w =0 

yyy . 




(79) 


for y s= 0 . Should these boundary conditions be given on other edges, 
the necessary changes in the formulas can be made easily. 

Several examples of the buckling of a simply compressed plate 
with various boundary conditions of the type just discussed are 
considered in the next section. In all these examples it is assumed 
that the plate is in the state of compression previously 
described, = -cr 0 , and that the edges x = a and x = 0, 

perpendicular to the direction of the compressive force, are simply 
supported. The other boundary conditions are specified for each 
example. Equations (75) are satisfied at x = Q and x = a if the 
deflection rate is written in the form 
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v ■=* f(y) sin (mnx/a) 


( 00 ) 


-where m is an integer* Thus a section of the plate in the buckled 
state, obtained by setting y « Constant, is described by a series 
of sine waves J the integer m gives the number of half waves. 
Substitution of this expression for v into the partial differential 
equation (72) yields the ordinary differential equation 


(aijt/aJ^-jjf - 2(mrt/a) 2 D 12 f il + T> 2 ^f ±r - a 0 h(m*/a) 2 f 


(81) 


in which the Roman numerals denote the corresponding derivatives with 
respect to y. With 


P 


2 



d 12 

®22 


2 ^ _ /gjA 2 ^ll^gg " S 12 g 

A -^12 V®/ DlpDpp 


y (82) 


equations (82) can be written in the form 


f iv - 2p2fii + p 2( p 2 _ „ o (83) 

The general solution of this equation in f(y) is 

f(y) * Qq G08h ry + c 2 sinil r y + °3 303 sy + 3 i), 3 ^ sy (84) 

where o-^, c^, c^, and c^ are arbitrary constants which must be 

determined from the boundary conditions, and 



NACA TN No . 1530 


35 


r = \j p(q + p) 
b = \Jp( q. - p) 


(85) 


Equation (8*0 is the fundamental relation -which must he studied 
for each particular case of buckling from a state of simple compression 
The boundary conditions of the type discussed lead to linear homo- 
geneous equations for the constants c-^, c q, c^, and c^ . The 


condition that these equations possess solutions c^, Og, c^, 

and c^ -which are not all zero yields an equation in r and s from 


■which the critical stress cr 0 can be determined. 


Specific Examples 

Plastic buckling of a narrow strlu; relation of present theory to 
beam theory .- It eoems -worth -while to investigate how the present 
theory of the plastic buckling of plates is related to the EngesBer- 
Von Karman theory of the plastic buckling of beams . Consider a 
rectangular plate of thickness h -which is simply supported 
along x = 0 and x = a, free along y = ±b/ 2 , and under the 
compressive stress a x = ~o Q . As b — > 0 for a fixed value of a, 

the buckling condition for this plate might be expected to approach 
that of a beam which has the length a, is simply supported at its 
two ends, and possesses a rectangular cross section of height h 
and width b. 

The solution given by equation ( 80 ) automatically fulfills the 
boundary conditions at x = 0 and x = a . In addition, the symmetry 
condition along y = 0 , equations (79) > requires that the coefficients 
in equation ( 84 ) satisfy 


rcg + sc^ = 0 


( 86 ) 


and 


r-^Cg - s^cj^ = 0 


(87) 
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Accordingly, 


c 2 = c h = 0 


(88) 


if the function f(y) is not to vanish identically. Boundary 
conditions, equations ( 77 ) > for a free edge along y = b /2 furnish 
the equations 


" £L 2-[pi2 " ( x " v )®]( c l COBh ^ + c 3 COB 

+ Dgg^Ci cosh - s 2 c 3 cos ^ = 0 

- + (1 ~ V)dJ (r Cl sinh ^ - sc 3 sin 

a 

+ siri11 + s ^ c 3 sin ~ 0 


(89) 


The boundary conditions corresponding to a free edge at y = -b/2 are 
fulfilled automatically if equations ( 89 ) are satisfied. Now these 
relations are linear homogeneous equations in c-j and c-^j if they 

are to have nonvanishing solutions, the determinant of the coefficients 
must be zero. Therefore, 




This det-enninantal equation can be reduced to the form 

(j 2 - m 2 k 2 E 2 q) q tan 3, + (q 2 + m 2 it 2 R 2 o) £ 


tanh ~ = 0 

e 


(9D 


where 


B = b/a 


i = rb s b\Jp(q + p) 
q = sb = b\/p(q - p) 

Q= [p^ - (1 - V)d]|d 22 


r (92) 


uo 
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If, for a fixed value of the span a, the -width b of the plate 
approaches zero, R and hence § and T) tend toward zero too* 
Accordingly, the functions 6 tanh (g/2) and T) tan (q/2) appearing 

in equation ( 91 ) may be replaced by | 2 /2 and t) 2 /2, respectively. 
For m = 1, in particular, this relation can be reduced to 


KDS2' “ * d 22' 2 ~ d 12' 2 + - 0 (93) 


in which the following notation has been used 



2 

« rr n 4 " 

k =. 

12 °l 

n 2 h 2 E ( 

Du 

' “ D 11 /D 

®12 ' 

’ - Dig/D 

D 2 2’ 

- D2 2 /D 

A = 

( D ll ' D 22 


(1 - V 2 ) 


(9«t) 


A detailed development of equation ( 93 ) is presented in appendix C* 
If the plate under consideration buckles within the elastic range. 
Dll' =» Dqg ’ = d 22 ' * consequently, A = 0 and Q = V. According 
to equation (93) then, k = 1 - V 2 and 


a 


o = 


JT 2 !^ 


12 a £ 


o 


(95) 
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The compress Its force P = a 0 bh under -which the plate buckles is thus 
seen to equal 


P 


n 2 h 3 bE o 

12a £ 



(96) 


where I ■ bh-Vl2 is the cross-sectional moment of inertia. 

Equation (96) is Euler's formula for the elastic buckling of a simply 
supported beam of span a. 

If the plate buckles after the compressive stress a 0 has 

exceeded the limit of proportionality, the evaluation of equation (93) 
becomes more difficult. Noting the relations given in equations (9k), 
equation (93) can be transformed into 


k 


D ll‘ 



(97) 


Eor v s this relation can. be handled quite sinqqly. A straight- 
forward computation (see appendix C) yields the following result for 

the critical compressive stress a = a : 

or o 


cr 


Ay i 

3a® (l 


The critical buckling load P is given by 


P 


C HV1 bh 

or 


^o 1 jfc 

a 2 (l + ix) 2 


(98) 


( 99 ) 


where again I = bh3/l2 is the cross-sectional moment of inertia. 
Hence the reduced modulus E* is given by 
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E* = 



(K>0) 


according to the definition of X given by equation (13). Tide 
reduced modulus is identical with that found in the Engesser-Yon Karman 
theory of the buckling of beams beyond the elastic limit. (See 
references 2 and. 3 0 Thus in the case of an incompressible material, 

V = — f the critical load for a beam can be found as a limiting case 


from the theory of plates . It should be noted that}- within the 
framework of beam theory, the critical stress is independent of the 
value of Poisson's ratio V. 


This result is not necessarily true for materials which are not 


incompressible. 



The quantity 


k _ 
1 - / 


(mi) 


which is simply a constant multiple of the critical stress for a given 
plate, has been evaluated as a function of X for the case Just 

mentioned,, v = and for a material with v = O.32; the results 

are plotted in figure 4. Although the two functions agree at X = 1, 
as previously proved, there is a marked difference which increases 
with increasing values of X. Consequently, two beams of the same 
shape but with different values of Poisson's ratio will buckle at 
different Critical stresses • 

At first sight, these results may seem somewhat contradictory 
since the analysis of plates developed herein is a generalization of 
the Engesser-Von Karmen theory of the buckling of beams and yet this 
theory does not always appear as a limiting case of the present 
analysis. An explanation of this discrepancy is afforded by a closer 
examination of the position of the neutral surface z 0 as a function 
of the parameter X • The function z Q depends on the plate 

thickness h and on V and X. The limiting process considered 
does not affect this relation and the neutral line in the beam is 
determined by the intersection of the neutral surface in the plate 
and the vertical plane y = Constant . The neutral line so determined 
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need not coincide -with the neutral line determined directly from the 
theory of beams. This difference is the basis of the discrepancies 

in k/(l - V^) noted. 

According to equations (i^), the function z Q is given by 


z 


o 


(2 - y) 6 1 4 - (2v - l)€ g 
K 


where 


K = (2 - V)^ + (2V - l)Kg 


The quantities and 6 g are the strain rates in the middle 

surface in the x- and y-directians, respectively. For the small 
displacements studied herein, is the rate of curvature of the 

line of intersection of the plane y = Constant and the middle 
surface . Similarly, Kg is the rate of curvature of the line of 

intersection of the plane x = Constant and the middle surface. In 
general, the position of the neutral surface depends on the four 

, , ♦ • 

quantities e^, K^, and K^ and this holds in the transition 

from plate to beam. On the other hand, the Inge sser -Ton Karman theory 
of beams assumes that the position of the neutral line depends on e-j_ 

and Kj_ alone. (See reference 1, p. 158 .) The position of the 

neutral line will depend on the type, of analysis used. Specifically, 
it can be shown (see appendix C) that the position of the neutral 
line found by considering the intersection of the neutral surface 
and a plane y = Constant and the position determined by the , 

Engesser-Yon Karman method agree if and only if V = i. Therefore, 

it is to be expected that the transition from plate to beam is valid 
only under this condition. 

Buckling of a simply supported plate .- The critical load for a 
simply supported plate under edge thrusts, stressed beyond the elastic 
limit, has been discussed by Timoshenko (reference 1, p. 3&7) on a 
purely formal basis. The object of this section is' the development 
of an analytical formula, for the critical stress of such a plate on 
the basis of equation (J 2 ) and the subsequent comparison with 
Timoshenko 's results . 
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According to equations (75 ) , the boundary conditions for a plate 
simply supported §t all four sidas can be -written in the form 


v = 0 

. 'PllV TT " Jj>22 " (3- " V )Dj Wyy =* 0 

at x s 0 and x = a, end 

v = 0 

” D 2g^yy ’ [ d 12 ' (1 - vjcjv^ = 0 


y (3X)2) 



h (103) 


at y = 0 and y = b* The boundary conditions at x *= 0 and I = a 
are fulfilled automatically if the function V is of the form given 
by equations (8 q) and ( 8 l) • Equations ( 103 )> applied to the 
function f(y) at y = 0, require that 


c-j_ *f- c = 0 

~%e(f 2c? i * + (jF) (pie ' ( 1 - v ^ D ](°i + c 3 ) - 0 


(iQfc) 


Consequently, 


°l = *°3 




0 


^ (105) 


thus 


o^_ = = Q 


(106) 


since r2 + ^ 0 if buckling is to take place. 
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The conditions at y = b Imply that 


C2 sihh rh +■ cj^ sin sb = 0 


-D22^ r ^ c 2 s thb. rh - s^C^ sin sh) 


> ( 107 ) 


[®12 ” D (1 - V)] (c2 sinh rh + sin sh) = 0 


which are equivalent to 


Cg slnh. rh + sin sh = 0 
r^c^ slnh rh - s^c^ sin sh = 0 


y (108) 


Wow these equations are linear, homogeneous equations in the two 
unknowns cv> and c^; they possess nontrivial solutions if and only 

if 


slnh rh 


sin sh 


r^ slnh rh -s^ sin sh 


= - (r® + s^) sin sh slnh rh = 


( 109 ) 


This condition will he satisfied if r^ + s^ = 0 , r = 0 , dr s’ = A'rf 
where n = 0, 1 , 2 , - . . • The first possibility is untenable; the 
second implies, according to equations (108) , that s = 0, n, 2k, . . . 
If s = 0 , then from equations (108) r must also vanish and buckling 
will not take place . The buckling condition must therefore be s = njt 
■there n = 1, 2, 3> * • • • 
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Since 

equations 


s = \/p (< 1 - p) 
( 82 ) that 


from equations (85), it follows from 


,2 _ S 2 ' .Z 

" b 2 
= p(q - p) 


mff \| D 12 \l a o h _ ^llPez ~ D lg 2 _ mir\pl 2 

a V I ^2 Y^lg V a ' ^ll 1 ^ a J^g 


(110) 


With the notation of the previous sections for Dqq’j Djjg Dgg 1 , 
and A, this relation can be written as 


n 2 ^ = mjr\ D lg 
b 2 a \ ®gg ' 


g o h _ /ffljA 2 P g2 ' A _ mjU D lg ' 
^DDig ' \ay D22 ' a ' D gg ' 


( 111 ) 


The critical stress 0 can then be found by direct computation to be 

cr 


°cr “ °o 


+ mb 

b 2 h \ m b V D 12' a 




p lg’ 

) + m 22 


D22' 

) 

Va/ 


( 112 ) 


Now each of the terms in the first parentheses on the right-hand side 
is positive; consequently, the minimum critical stress will occur 
when n = 1 . Thus the critical stress is given by 
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a 


cr 


it^DDjg' /i a\l I> 22 ' 

A ^ b Pl2’ 




(113) 


In the elastic case, X = 1 and the functions appearing in the 
formula for the critical stress take on the values 1 = P~)p ' = Dgg ’ ~ 

A = 0; consequently, 

(lUO 



This result is the same as that found for the elastic case directly 
(reference 1, p. 330 ). Timoshenko suggests (reference 1, p. 3®7) the 
application of the following equation for the "buckling of a simply 
supported plate compressed "beyond the elastic limit. 


i (aw 


aw J u J + 
im 


2 \fa w. 


xxyy 


+ w_ 


yyyy 




N w 

X XX 


= 0 


(115) 


where here a = IS*/. E Q , the ratio of the reduced modulus (defined in 

equation (IDO)) to Young's modulus. Under these assumptions, the 
critical compressive stress is found to "be 


a 


or 



(116) 


Although Timoshenko (reference 1, p. 3^7 ) gives only the critical 
stress for the half “wave number m = 1, his results can be easily 
extended to the form given in equation (ll6) . 

It will be found more convenient to use a new quantity k rather 
- hVip-n <y for comparison of the two theories . This parameter is a 

modified form of k defined in equation (9*0; more precisely. 
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k6 


'-S* 


2 

12a b 
cr 

it 2 ! 2 !. 


C 1 - v£ ) 






(117) 


For the theory proposed in the present report, this quantity becomes 


k 


Bj2 



I> 22 t An .% 2 


(118) 


■where E = b/a. 


as before- In the elastic case, TE takes the form 


k = 



+ mR 




(119) 


On the other hand, Timoshenko 's result can be written as 


k = 



( 120 ) 


The differences in the two results appear more readily if 
equations (ll8) and (120 ) are expanded in full. After same simplifi- 
cations, equation (ll8) becomes 


k = 


D u 'i!i¥ + 2Dj2 


D 22 ' 

m% 2 


( 121 ) 
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■while equation (120) takes the form 


k = + 2 \fa + — 

m % 2 


( 122 ) 


It is seen then that k is the same type of rational function 
of for "both theories except for the fact that the coefficients 

are different functions of X - 

Unlike in the case of the cruciform section which is treated in 
the next section, the confutations involved in evaluating k as a 
function of X and mR are relatively simple since the coeffi- - 
cients 1 , D ^ ' , and a can he tabulated once for all 

independently of the stress -strain lav and the particular geometrical 
ratios under consideration- The critical values of Tc can be found 
with reasonable speed by the following procedure. Curves of k 
against X can be obtained for various values of the parameters m 
and E by evaluating equation (121) or (122) , depending on which 
theory is used- These results, which form, the bulk of the computation, 
do not depend on the stress-strain law but only on V. Consequently, 
the curves of TE against X thus found are valid for_all materials 
having the same Poisson's ratio- On the other hand, k and X are 
related by a second equation which depends on the stress -strain law; 
namely, from equation (117), 


k = 




(1 - ■£) 


Since a is a function of X, determined by the given compressive 

stress-strain law of the material in question, k also can be plotted 
as a function of X once the plate width-to-thickness ratio b/h 
has been fixed for a particular example - If these curves are plotted 
on the same sheet as those previously described, the intersections 
will give the critical value of k, and hence the critical stress, 
corresponding to a given stress-strain law and_given plate param- 
eters b/h and E- It should be noted that k depends on the 
ratio b/h in the plastic range, whereas this is not the case within 
the ordinary elastic theory. 
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Since the quantity X enters into the confutations for the 
second set of curves through the stress-strain relations, it is quite 
useful to represent the stress-strain curve analytically . An 
expression of the form. . 


6 


==- a + Ac 


(223) 


has been fitted to the experimental data in such a way that the 
experimental and theoretical curves pass through the same initial 
and final points and the slopes at these points coincide. The fitted 
curve together with the experimental points is shown in figure 5 • 

The computational program described has been carried out for the 
material given in figure 5 and the ratio of plate width to plate 

thickness fixed at b 2 /h 2 = 1000 • The solid curves of figure 6 
give k as a function of X with parameter B as defined by 
equation (121), each set of curves representing a different value 
of m. These curves can be applied to a plate with any ratio of b/h 
and any compressive stress-strain law provided, v = 0- 32. The dashed 
lines, on the other hand, represent the relation between TE and X 

given by equation (117) > they correspond to b 2 /h 2 = 1000 and the 
stress-strain curve of figure 5* The points of intersection have 
been determined and the critical values of k plotted in the usual 
fashion as functions of 1/te in figure 7* Timoshenko's results are 
also shown as well as those for the elastic theory. Since the critical 
stress is given by TE multiplied by a constant, the interpretation of 
these curves of TE against l/R holds equally well for the critical 
stress as a function of l/R . 

First, it should be noted that the plate will buckle in one or 
more half waves according to the magnitude of the ratio l/R = a/b; 
this holds for plastic buckling as well as elastic buckling. 
Furthermore, it is necessary only to consider that part of the curve 
corresponding to a given value of m which lies below the inter- 
sections with curves belonging to adjacent values of m. Under these 
conditions, it can be readily seen that the critical stress obtained 
from the present theory lies between the results of the elastic 
analysis and Timoshenko's formal procedure. More precisely, the 
elastic critical stress is higher than that predicted by the present 
theory, whereas Timoshenko's buckling stress is lower- The transition 
from a buckling mode with a given number of half waves to the next 
higher occurs at practically the same values of l/R in both plastic 
theories ) whereas, this transition occurs at slightly larger values 
of l/R in the elastic range • Finally, the locus of m i n i m ums of 
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all three sets of curves Is, to a high degree of approximation, a 
straight line. While this straight line is fixed for all values 
of b/h in the elastic case, it will shift in the plastic case as 
this parameter is changed. 

With a slightly different computational procedure, the results 
Just obtained for the buckling of a simply supported plate by the 
present theory can be compared with Ilyushin 's solution of the same 
problem (reference 6) . Ilyushin 's general equation for the buckling 
of a rect ang ular plate compressed in one direction is of the same form 
as equation (72) except for the difference in the coefficients . His 
fundamental relation can be written as 


J 11 xxxx 




(1210 


where the coefficients D^*, D^p*, and p * are rather compli- 

cated functions, of the stresses which will be discussed shortly . 

Since these results apply only to incompressible materials, the theory 
developed in the present report must be specialized to the case for 

which V = i. 

For the case of a simply compressed plate, the basic functions 
entering in the definition of the coefficients in equation ( 1210 can 
be given in a slightly less complicated manner than that used by 
Ilyushin for the more general problem. Let E s denote the "secant 

modulus, " that is, quotient of the compressive Btress divided by. the 
compressive strain as obtained from a compression test for the 
material in question. Then Ilyushin's function co can be shown 
to be 



(125) 


where E Q is Young's modulus for the material- With Von Kaman's 
modulus, E*, as given in equation (100), a new quantity ic may be 
defined by setting 


k = E*/E 0 


(126) 
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The function •will "be used to designate the following relation 
containing k and 


* = 




(127) 


The coefficients D^*, and Dgg* can then he shown to he 

of the form 


D u* - ^ - ♦’C 1 - i - j - ? ir) 

D^* = D(1 - t) 


y (128) 


d 22 * = D(1 - *) 


where D is the flexural rigidity of the plate defined in 
equations (73) • It should he noted here that Ilyushin's coefficients 
depend on the two moduli E and E s , whereas those appearing in the 

present theory depend only on the tangent modulus E or rather the 
ratio X = E/E 0 . Consequently, it is extremely difficult to carry 

out the larger part of the computations for k independently of 
the stress -strain law as can he done with the method presented herein- 
The numerical technique must therefore he changed somewhat- 

The boundary conditions for the simply supported plate in 
Ilyushin 'e analysis are 


w = 0 


- D ll**xx - i°ss **yy = 0 


at x = 0 and x = a, and 


(129) 
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w = 0 

“ rp2t?^xx ~ D 22*^yy = 0 


(130) 


at y = 0 and y = b. These boundary conditions are completely 
analogous to equations (102) and (203), provided Poisson's ratio V 
is taken to he l/2_in the equations (102) and (103)* If the new 
quantities D.^, D ip , and I) 2g are introduced in the following 

manner 


P 11 “ D U*/ D 
d 12 = D ] 0*/ D 
^22 = D 22*/ D 


(131) 


the critical parameter Tc can he found in precisely the same fashion 
as that used in developing equation (121). Consequently, 


k = D 11 m £ R 2 + 


2D 12 + 


22 


m. 


L% 2 


(132) 


As has been previously noted, the coefficients 35^, D-^, 

and Dpp are functions of both the secant modulus and the tangent 

modiolus with the result that specific reference to the stress-strain 
law in compression must be made in order to evaluate equation (132) 
readily. The critical quantity TE has been computed as a function 
of cr™ for the material shown in figure 5 with m and E appearing 

vx __ 

as par ame ters ■ For a given width-to-thickness ratio b/h, k can be 
determined by finding the intersections of these curves with the 
straight line defined by equation (117). This procedure has been 
applied to the case b 2 /h^ = 1000 previously discussed. The results 
obtained from Ilyushin's method and from the present method 

with V = ^ are shown in figure 8. Again considering only those 
parts of the curves which lie below the intersection points corresponding 
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to consecutive values of m, certain general conclusions can be 
drawn. The minimum value of k as determined by Ilyushin's method 
is practically constant for all values of m as has been seen in the 
other theories considered. On the other hand, Ilyushin's theory 
predicts smaller values of k and hence lower critical stresses 
than those of the present theory. Comparison with figure 7 would 
indicate that Ilyushin's .critical stresses lie roughly between those 
of Timoshenko and the present theory. Finally, buckling occurs at 
practically the same wave number in both theories although the jump 
in form takes place at slightly smaller values of l/R in the 
present theory . 

Torsional buckling of cruciform sections .- The torsional buckling 
of a cruciform’ section under compression can be studied by treating 
each flange as a simply compressed plate, simply supported at y = 0 
and free at y = b. (See reference 1, p. 340.) As in the previous 
case, the deflection rate- w will be given by equation (8o) where 
f (y) Is found from equation (81) . The conditions for the simply 
supported edge (equations (75)) at y = 0 become 


c 3 = ' C 1 
f i:L (0) - 0 


(133) 


The last equation implies r^c^ = s^cg. Since c^ = -C]_ and r 2 
and s 2 do not vanish, cj_ = c^ = 0 ; consequently. 


f(y) = °2 slri11 + c j| sln B 7 


(134) 


Finally, the boundary conditions at the free edge y = b require, 
according to equations (77 ), that 



- (l - V)dJ(c 2 sinh rb + c^ sin sb) 

+ sinl1 “ s2 °4 sin s ^) = 


0 





+ (l - V)dJ(tc 2 cosh rb + sc^ cos sb) 
+ ^ 22^^ c 2 co ®k r 4 " s ^ c 4 cos = 0 


( 135 ) 
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Nov equations (135) are linear homogeneous equations in the two 
unknowns eg and c^. If these equations are to yield nonvanishing 

solutions, the determinant of the coefficients must vanish; tha t is. 


a l b l 


a g h 2 


= 0 


(136) 


where 


a l = ^22 r2 " ^ m 2 lt 2 §i2 - 0 - " v) 5 | /a 2 J^ 

b l = ‘ 2® 2 + - C 1 " v )d] /el 2 

a 2 = r^Dggr 2 - -jm 2 * 2 ^ + (l - V)dJJ a 2 ^ J 

h 0 = -b^D 22 s 2 + -|m 2 n 2 [^22 + (l - v )®] / & 2 J^ 


sinh rh 


sin sh 


cosh rh 


cos sh 


a^tematic interchange of the hyperbolic and trigonometric 
functions, the determinant in equation ( 136 ) is the same as that 

appearing in equation (90). The determinantal equation can then 
he reduced to 


W* - - a - *>§ 


tanh rh 


= r -< 


Dggs 2 + sE^jp-^ - (1 - V)d] 


r tan sh (137) 


With the notation introduced in equations (92), this result 
rewritten in the form 


can he 
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(i 2 - m 2 A 2 «) 2 tofi. . (* + ,2,%^ Jas-a (138 ) 

Again I and q can "be determined as functions of the quantities Du', 
D ip ', Dgg and k defined in equations (9*0* With some additional 

transformations, equation ( 13 ®) can then "be put into a form suitable 
for solving for the quantity k as a function of X • Although the 
computations are somewhat tedious, the solutions obtained are completely 
independent of the stress-strain law. These results are represented by 
the solid lines of . figure 9 and are marked for various values of the 
side ratio R. The dashed curves, on the other hand, represent the 
quantity k as a function of X as derived from equation (117)* 

These curves have been computed for the stress -strain law shown in 
figure 5 and the values of the vidth-to-thickness ratio b/h as 
indicated in the graph. Given the two dimensionless quantities b/a 
and b/h, the critical stress is then determined by the corresponding 
solid curve and dashed curve . This point giveB the desired value 
of k, from which the buckling stress is found by solving 
equation (117) • The details of the procedure are found in appendix C. 

It should also be pointed out that the solid curves have been computed 
for the half-wave number m = 1, for it is shown in appendix C that 
the lowest value of k, and consequently the lowest critical stress, 
will be obtained for this value of the wave number. 

Certain general conclusions concerning the buckling of such a 
section can be reached without reference to the stress-strain law by 
means of the solid curves of figure 9* These curves show that, for 
a given side ratio R, IE is a decreasing function of X • This 
decrease is not the same for all values of the side ratio, however. 

For small values of R (large length- to -width ratios), IE is almost 

constant. In fact, when R ^ 0.10, k is practically constant and 
the critical stress is the same as that found in the theory of elastic 
buckling of plates . The plastic effects are very pronounced fur 
larger values of R and increase as R — > 1, that is, as the 
rectangular plate becomes more nearly square. 

Finally, the intersections of the solid and dashed curves of 
figure 9 yield k as a function of the side ratio R; these are 
replotted in figure 10. The top curve represents the elastic case 
(X = l) and is independent of b/h. The other three curves 
correspond to the cases b^/h^ = 300 , 250 , and 200, respectively. 

It should be noted that k for the elastic case is always greater 
than k for the plastic case; this difference is more marked for 
shorter (larger R) and thicker (smaller b^/hS) plates. 
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Energy Method . 

In the preceding section rigorous solutions were obtained for 
several cases of plastic buckling of simply compressed plates. In 
cases for which the rigorous solution becomes too unwieldy or is not 
known at all, the energy method provides a convenient means of 
finding approximate values of the buckling loads . 

As far as the elastic buckling of plates is concerned, the energy- 
method is well established. (See, for instance, reference 1, 

PP- 325 ff -) At first glance, it may appear somewhat doubtful 
whether this method can be extended to the plastic buckling of plateB. 
That such an extension is legitimate, however, follows from the fact 
that equation ( 109 ) has the same form as the equation for the elastic 
buckling of an orthotropic plate with flexural rigidities D-, n 

and D pp and torsional rigidity D.^. Since the energy method is 

valid for this problem of elastic buckling and since the plastic 
buckling of a simply compressed plate is governed by the same equation, 
the application of the energy method to the plastic buckling of plateB 
Is legitimate. 

Under sufficiently small edge thrusts the flat form of the plate, 
represents a stable equilibrium configuration. When the edge thrusts 
reach the critical value, however, this equilibrium becomes indifferent 
and the plate may assume a bent form. This transition from the flat 
to the bent form, that is, from one indifferent equilibrium configura- 
tion to another one, does not involve any energy input or output . 
Accordingly, the work done by the edge thrusts equals the flexural 
energy of the bent plate • The mathematical expression of this 
principle is the basic energy equation: 



D f j J [l - cS(2 - V) 2 ] w^ 2 + 2[v - cS(2 - V) (2 V - l)] w^v^ 


+ [1 - C»( 2 V _ + 2(1 - 


dx dy 


c 0 h j j w T c dx dy 


(139) 


where the range of integration Is the area of the plate, that is, 

0 ^ x a a, 0 = y = b. As is easily seen from equations (^6) 
and (7*0, the left-hand side of equation ( 139 ) equals 
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■which represents twice the rate at which work is done by the bending 
and twisting moments- In the elastic range, c = 0, this expression 
can be written in the familiar form 


(*xx + *yy) 2 “ 2(1 " v) (wyy " ^xy 2 )J *= & (l4l) 


The right-hand side of equation (139) represents twice the rate at 
which work is done by the edge thrusts and has the same form as in 
the elastic case. 

The following expression for the critical load of the plate is 
obtained from equation (139) : 



*o h = s' 


(142) 


where 



A = D II ^ C 1 “ c5(2 “ V)2 ]"xx 2 + 2 [ v ~ °5(2 — V ) ( 2V - 1)] w w 


"W 

XX yy 


+ [l - c5(2V - l)] w 2 + 2(1 - v)w 2 > dx dy 


(143) 


and 


-//t 


B = / / v x dx dy 


(144) 


The right-hand side of equation (142) depends on the function w, 
that is, on the deflected shape of the plate. Now, any restriction 
which is imposed on the deflection rate w, over and above the 
boundary conditions discussed in the third section of the ANALYSIS 
amounts to an increase in the stiffness of the plate and must, 
therefore, lead to an increase in the critical load. The critical 
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load for a given plate is accordingly found as the smal l est value 
which, the right-hand side of equation (14-2) can assume for functions w 
possessing continuous partial derivatives of the second order and 
satisfying the boundary conditions for the deflection of the plate . 

In the elastic case, c = 0, the right-hand side of equation (l42) 
is independent of the buckling stress a Q • Equation (142) therefore 

furnishes the buckling stress a D which corresponds to a given plate 

thickness h. In the plastic case, however, the numerator A of the 
ri^it-hand side of equation (l42) depends on the buckling Btress cr 0 

through the quantities c and 8. Equation (142) should therefore 
be considered as an equation which determines the critical thickness h 
corresponding to a given buckling stress . If the expression for D, 
from equations (73) , is substituted into equation (l43) and the 
resulting form of equation (142) solved with respect to h, the 
following formula is obtained: 




12cro (l - v2)b 




(145) 


where B is given by equation (144) and A 1 denoteB the integral on 
the right-hand side of equation (l43) • 

T-he energy method will now be applied to the buckling of a 
cruciform section previously considered. A suitable choice of the 
approximate deflection rate w is 

V = q(i - 2| 3 + I 4 ) (146) 


where | = x/a and q = y/b. This function satisfies all the boundary 
conditions on the plate except at the free edge y = b • Nevertheless, 
this expression proves to be a satisfactory appro xima tion as will be 
seen_shortly • The necessary integrations indicated in the definitions 
of A and B (equations (143) and (l44)) can be easily carried out 
in this case; for example. 


J U a do Jo 


1.6R 

a 2 


(147) 
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dx dy 




C 1 


12t 2 + 8a 3 + 365^ - 48| 5 + 16 a 6 ) dl dri 


17 

353%? 


(148) 



dx dy = B 




- 12I 2 + 8a 3 + 36I 4 - 48I 5 + i6S^) dl drj 


- 17B 
105 


(149) 


where B = b/a. Thus 


_B 

A' 


UepR 2 - 


i68r‘ 


Cl- 


c5(2 - V) 


2 ]* 


102(1 - V) 


17a % 2 


168*?%^' + 102(1 - V) 


(150) 


and 


k = 


le^D^' + 102(1 - V) 
17* 2 


(15D 


where k is defined in equation ( 117 )* 

Since D-^' 4as 4een previously computed as a function of X 

(see table I), k may be readily evaluated from equation ( 151 ) as a 
function of X for various values of the parameter B . This has 
been done for the example V = 0 >32 considered in an earlier section. 
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The results are shown In figure 11 in ■which the solid lines represent 
the solution obtained by the energy method and the dashed lines are 
those found from the exact solution. A brief inspection will show 
that the energy method, as applied here, gives a very good approxi- 
mation in the technically interesting region of small values of E • 

The error increases as E becomes large and reaches a maximum of 
about 8 percent. It should also be noted that for a given value 
of E the error is an increasing function of The energy method 

actually yields better results when the complete problem of deter- 
mining the critical values of k associated with a given stress - 
strain law is carried out. These points are found from the intersections 
of the curves of figure 11 and the curves derived from the stress- 
strain law through, equation (117) (fig* 9) • Since the curves of 
figure 9 are monotonic, increasing functions of X with a slope 
angle of less than 90 ° , "the error will be smaller than the original 
estimate. The results of applying this method to the cases treated 
in the previous section on exact analysis are shown in figure 12. The 
solid curves represent the solution by the energy method, whereas the 
dashed curves are the exact solution. It is seen that the maxi -mum 


error is of the order of 2— percent of the value of 

2 


k. This appears 


to be well within the order of accuracy of the theory itself. Finally, 
the approximate value of k is always greater than the exact solution 
as is usually found in the application of the energy method. 


Equation (151) also affords an easy means of studying the general 
character of the critical parameter k as a function of X for 
various values of E- It has already been pointed out that ' is 

a decreasing function of X starting at D-^' = 1 for 1=1 and 

approaching Du’ = 0 •2^13 for X — > co. Thus for small values 

of E, the first term has little influence and TE is approximately 
a horizontal line defined by the equation 


% 10g(l - V) 

17 * 2 


(152) 


The elastic solution based on the energy method .with the deflection 
given by equation ( 146) is 


£ = i 68 r 2 + 102(1 - v) 


(153) 
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and equations (152) and (153) agree quite closely for small values 
of E. Under these circumstances, then, the critical stress deter- 
mined from the plastic theory will he the same as that obtained 
from the elastic formulas. 


Brown University 

Providence, R. I., March 26, 19^7 
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APPENDIX A 

DIEFICUIillES PRESENTED BI ILYUSHIN'S STEESS- 
SIEAIN RELATIONS 


It has teen pointed out in the first section of the ANALYSIS that 
certain objectionable discontinuities are present at the neutral 
surface -when the state of stress and strain in a body is described by- 
relations of the type used by Ilyushin (reference 6) . More precisely, 
it has been shown that the equations for loading and -unloading do not 
give the same increments in strain at the neutral surface; this 
difficulty is not present in the theory of plastic flow developed 
in the present report. Such problems always arise when an attempt 
is made to take unloading into account in a theory of plastic 
deformation. (See reference 10, p. UOO.) 

The inconsistency is brought out more clearly if the whole 
problem is considered in terms of the state of stress and strain 
existing in the buckled plate- Let 


0 



co 

1 - CO 


(Al) 


Then Ilyushin's stress-strain relations for loading in the case under 
consideration (equations ( 30 )) can be written in the form 


2E 0 e x = (1 + fi)(2cr x - cr y ) 

2E 0 6y = (14-fi) (-q x + 2d y ) 

Vxy " 3d + 0)*^ 

Differentiated with respect to time, these become equations of the 
following type: 


(A2) 


“o 4 ! - d + n >(e*x - °j) + - Oj) 


(A3) 
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In particular, for buckling from a state of simple compression, 
c x = -cr 0 and cr y = = o, 


fi = o ± 
cLc?! 


- -f n '( 2 ^ - “y) 

where Q' = dQ/da^- Accordingly, 

>*x = (l + + *o n ) (2d x - «J y ) 


(Ak) 


2E, 


(A5) 


The boundary between the regions of loading and unloading is a 
plane z = Constant. If the displacements are to be continuous across 
this plane, e x , £ , and y , too, must be continuous. Moreover, 

according to Ilyushin's theory the sign of -o Q £ x serves as a 
criterion for loading and unloading. Since e x is continuous, it 
must vanish at the boundary between the regions of loading and 
unloading. Thus, 2d x - c y is continuous, too, and ft = 0 on the 
boundary. Denoting the values on the loading side of the boundary 
by the superscript + and the values on the other side by the 
superscript -, the strains at the neutral 

2E 0 e x = (1 + fl)(2a x + 

= 2 c x - CTy 

= 0 

2B 0 6y = (1 + «)(-a x + 

- + 2 V 

e 7 —— = 3(i + + 

o sy s xy 


surface become 


- V) 


+ 2cr. 


V) 


(a6) 
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Thus 


cr^ = (l + ft)cr I + 

V - (1 + n) V 



(1 + n)+ 




(A7) 


Since (l + fl) > 1, the rates of stress on the unloading side of the 
boundary are of the same sign as those on the loading, furthermore, 
the absolute values of the rates of stress on the unloading side are 
greater than those on the loading side . While this prediction of 
Ilyushin's stress-strain relations does not violate the equations of 
equilibrium., it does seem, quite strange -when compared -with the usual 
notions of loading and unloading. The result is of a sufficiently 
startling character to call for direct experimental verification 
before a theory of structural stability in the plastic range iB based 
on these stress-strain relations. 

That the discontinuities implied in Ilyushin's stress -strain 
relations are a source of concern to Eussian scientists working in 
this field, is seen from the following passage which is quoted from 
W. W . Sokolovsky 's recent paper (reference 14) in which the stability 
problem Is specifically discussed: 

"In solving these problems, the Mises-Hencky theory was employed, 
under the assumption of the incompressibility of the plastic material. 
This assumption has a strong influence upon the value of displacements 
and does not permit satisfying completely the conditions of the 
continuity of all components of stress and displacements on the 
boundary between the elastic and the -plastic zones ■ An objective of 
recent work is to avoid these defects of the theory." 

While the discontinuities appear to have been noticed, the 
source has not been traced correctly, for theBe discontinuities 
subsist even if the compressibility of the material is taken into 
account. They can be avoided only by replacing Ilyushin's theory of 
plastic deformation by a theory of plastic flow of the type given in 
the present report . 
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APPENDIX B 

DETAILS IN DEVELOPMENT OE FUNDAMENTAL PLATE EQUATION 
Determination of Neutral Surface 


The neutral surface can he determined from the condition that the 
additional stress resultants present in the buckled state must be in 
equilibrium. Rather than confute the rates of the stress resultants 
arising from c x and dy separately, it is much simpler to consider 

the combinations a x - Vdy and dy - Vo^. (A similar device is 

employed in reference 6 , p. 3^1*) For unloading, equations ( h- 0 ) 
may be used to yield 


r (bd 


On the other hand, equations (39) give for loading 



€3. = a x - vdy + (\ - l)cr z “ X 2 " 

• •. X ~ 1 • X - 1 V 

6 y - °y - Vct x " —2~ a x + “IT" a y 


(B2) 


Transposition in the first of equations (B2) leads to 


°x - Vc V - «x - ( 2b x ' *7) 


(B3) 


■which, according to equations (hi), can be written as 


O’ — V O’ = € 

3C y X ( er _ _ ( ~\ _ 




(5 - 4V)\ - (1 - 2V) 
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With 


c 


X - 1 

(5 - bv)\ - (1 - 2 v ) 2 


(B5) 


and K and z Q 


as defined in equations (h9), equation (Bh) becomes 


a x - V 0 y = e x + 2 cK(z - z 0 ) 


(B 6 ) 


Tlie quantity c depends only on the properties of the material 
through V and on the compressive stress through X • It is a 
nondecreasing function of X, for 


dc = (5 - hvIX - fl - 2V) g - (5 - hVHX - l) 
dX ‘ £5 - 4v)X - (1 - 2V) 2 ] 2 

= 5 - 4v - 1 + hv - hv 2 

[(5 - ^v)X - (1 - 2V) 2 ] 2 

_ b(l -v 2 ) 

[(5 - bv)\ - (l - 2 v)^J 2 

= 0 (B7) 


since -1 = v = i. For X = 1 (the smallest value of X), c = 0; 

for X = co, c = l/(5 ~ hv) . A similar espression for ay - Vc x can 
be found by the same technique; namely. 


- VCT X = 


€y - CK(z 


(B 8 ) 
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0 * 

The rates of change of the stress resultants N x and N y can he 

computed from satiations (Bl) , (B6) , and (b 8) . Two cases arise 
depending on -whether the region z,> z Q is a region of loading or 

■whether z < z Q is the region of loading, that is, -whether K > 0 

or -whether K < 0 . 

Case (l), K > 0 .- Since K > 0, the region z > z Q is the 
region of loading and z < z Q that of unloading. Consequently, 

equations (Bl) are -valid for z < z Q and equations (B6) and (b8) 

> ' 
hold for z = z Q . Now the rates of the stress resultants N x 

and Ny are defined as 



According to equations (Bl) and (B6)‘then 


(B9) 



(BIO) 


The strain rate e x = - zK from equations (b5) ; neither the 

function e ^ nor K depends on z . Therefore 



dz 



e l h 


(Bll) 



NACA TN No . 1530 


6 ? 


and 


- VNy = €* 


h + c. 


2 


In the same way, it can. be shown that 


(B12) 


N, 



■y 




(B13) 


From, equations (Bll) and (B 13 ) it has been shown in the second section 
of the ANALYSIS that the position of the middle surface is given by 
that value of z Q which is the solution of equation (58) • "With 

f; o = 2z o /h, equation (58) becomes equation (60) ; and with equation (62) 
the solution can be written as 


£0 = a * 




(B 14 ) 


Since c is a nondecreasing function of X with values in the 
range 0 = c = 1/(5 “ 4 v), the quantity a will lie in the range 

-co < 0 < -1 (B 15 ) 


Consequently, only the positive root in equation (Bit) can be kept 

if the inequality J^ Q | = j 2 z Q /hj =1 is to be preserved. It will 

be convenient to denote this root by £ + : that is, 

o 

£ o + = = a + \Ja 2 - 1 (Bl6) 

This is the formula for the neutral surface given in the second 
section of the ANALYSIS. 
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Case (2) , K < 0.- The same technique, as used for case (l), with 
the exception of certain changes in sign can be employed in the case 
when K < 0. "When jt < 0, the region z < z Q is the region of 
loading and z > z Q that of unloading. Consequently, equation (Bl) 
is now valid for z ^ z Q while equations (b 6) and (B8) hold 
for z = z Q . Thus N x - VN y becomes 

N x - VN y = P h//2 e x dz + 2cK P Z ° (z - zj) dz (B17) 

J-h/2 J-h/2 

which can be evaluated, as in the previous case, in the form 


K - vir y “ *l h - <*(-§ " z o) 2 

= Sjh - cK^| + z D ) 2 (Bl8) 

Similarly, 

^y " v3 *x “ *2 h + |^(| + z o) 2 (B1 9) 

• • 

Again the rates of the reduced stress resultants N x and N y must 
vanish. This gives 



(B20) 


Equations (49) can be used to show that 
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z c = (2 - V)i + (2V 
IL 

- - ^ (5 - *) 

K 

and equations (B20) and (B2l) yield 




0(5 - Iv) 


(B21) 


(B22) 


Equation (B22) is the analogue of equation < 58 ) which was found for 
the case K > 0. It can also tie solved by introducing the nondimen - 
sional quantity £ 0 defined as = 2^/h. The quadratic equation 


for is 


(1 ♦ E 0 ) 2 - 




c (5 - lv) 


= 0 


(B23) 


and this can he transformed into 


So +* 


£ ' c(5 - 4v)] 


+ 1 = 0 


(B2l) 


The quantity appearing in the brackets was defined as a in 
equations (61 ) } thus the solution of equation (B2l) is given by 


t = -a ±\ja? 


(B25) 


Again, the restriction that |^ Q | =1 requires that only one of the 

signs in equation (B25) can be taken, namely, the ml nu s sign- Thus 
the solution is that given by equation ( 63 ); namely, 


£0 = 


-a 


- \ja? - 1 = £ c 


(B26) 
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where the symbol £ 0 has been introduced to distinguish this case 
from the root £ q + previously found. Comparison of equations (Bl6) 
and (B2 6 ) shows that 



(B27) 


In addition, £ 0 + < 0 for E > 0 and £ 0 >0 for E < 0, as 
might be expected from the geometry of the situation. 


Determination of Rates of Change of Bending 
and Twisting Moments 

The moment rates can be computed now that z Q or is known 

a • 

for each of the cases E > 0 and E < 0 . Since the rates of change 
of the reduced bending moments and My are defined as in 

equation (64), the device of computing M^. - VNL. and My - V M X may 
be used here. 

Case (l) , E > 0 .- Application of equations (Bl) and (B6) yields 



The quantity z q 


can be replaced by 



with the result 


Mjj. - VMy = 



<=£(1 - §5 0 + 



(B30) 
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With 


28 = 1 


-| ? 0 


+ + 1 c +3 
2 ° 


the last equation becomes 

- v «y - - (% " 2c8 ^) 


Equation (b 8) for -Vcr x differs from, equation (b 6) for a x 
only in that -K replaces 2K and £ y replaces £ x in 
equation (b 6) . Without further computation, it can be seen that 


My - Vli X + C5i ) 


Equations (B32) and (B33) can be solved for the reduced bending 
moments -with the result that 



(B31) 


(B32) 


(B33) 


(B3^) 


(B35) 


The curvdture K can be expressed in terms of K^ and Kg according 
to equations (i<-9) ; consequently. 



M I =- 


12(1 - V 2 ) 


K X [l - c8(2 -V) 2 ] + K^[v - c8(2 - V)( 2 V - 1)] 


(B36) ^ 


J 


12 




Kpjv - c8(2 - V)(2V - 1)] + E^jl - c8(2V - l)^j 


(B37) 


These relations are the same as equations (67) and (66) -which were given in the second section 
of the AWAIZSIS without proof. 

Case (2)^ K < 0 .- This case can he handled in exactly the same manner as the previous one 
provided the integrals are split up as follows: 


*r- 



e z 
x 


dz + 2cE 




o 

h/2 


( z “ z q) 


dz 



« 

For E<0, 



and 




K, 



(B 38 ) 


(B39) 


(Bho) 
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No-w 5 Q ~ = - 5 0 + ; thus, 

( 2 . - 1 5 " 3 + 1 Z ^ = (1 + I 5 +3 -15 +N ) = 25 (BUl) 

V7 2 o 2 o ) \ 2 O 2oy 

In other words, M x - VMy will he given hy the same formula for K < 0 
as for K > 0. Similarly My - VM X will he given hy the same 
expression in either case; and, consequently, and My can he 

found from equations (B 36 ) and (B37) independently of -whether K < 0 
or K > 0 . 

The rate of change of the reduced twisting moment M^y is 

defined in equation ( 65 ) • According to equations ( 40 ) and (4-5) 
equation ( 65 ) may he rewritten as 


Ph/2 . 

V - - a i ri - v T , 42 

U -h/2 


ay 


(B42) 




(B43) 


Integration of the last equation furnishes 


■ _ h 3 (1 - v) £ 

^ -(1 - v 2 ) 2 12 


12 


(B44) 


This result, too, holds independently of whether K > 0 or K < 0 . 
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APPENDIX C 

DETAILED COMPUTATIONS FOR SPECIFIC EXAMPLES 
Buckling of a Thin Strip 


According to the definition's of r and s, equations (85), 
and of p, equations (82), 


* . S 2 . 2p 2 = La 


W D, 


(Cl) 


'22 


The terms of the second row of the determinant (equation (90)) can 
therefore he written in the form 


j^s 2 - *^$12 - (1 - V)D]| slnh f 


®22 r2 - ^[fl 2 ~ C 1 " V ) D ] 


sin 


Bh 


( 02 ) 


The determinantal equation (90) is thus seen to he equivalent to 


-Td^ . - (X - V)d]J- tan * 

+ ^£[ Dl2 - (1 - v ) n ] j . 


+ r 


tanh ~ = 0 (C3) 
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With, the definitions of equations (92), equation (C3) can he 
re-written in the form of equation (91) ■ The evaluation of the roots 
of the transcendental equation ( 91 ) can he simplified further by 
introducing equations (9*0 • A simple computation -will show that 



(C» 


If, for a fixed value of the span a, the width h of the plate 
approaches zero, R and hence £ and q tend toward zero too. 
Accordingly, the functions £ tanh (t /2) and q tan (q/2) appearing 
in equations (Ch) may he replaced hy £ 2 /2 and q 2 /2, respectively. 
For m = 1, in particular, equation ( 91 ) can then he written aB 
follows : 


kD 22 ' " AD 22 ,2 - D22 ' 


') 


(\|lcD 22 ' - AD 22 ' 2 + ' - 01*22^) (\[ 

+ QmD 2 2 7 - AD 22 ' 2 - Di2 ' + QD22^) 2 (^22 ' “ ^°22 ' 2 + D 12 ^ = 0 ( c 5) 


or in the form of equation (93) as 


H )22 ' - ^22 ' 2 " d 12 ' 2 


+ Q 2 1 ) 22 ,2 = 0 


Equation (99) for the buckling load when V = — can he 

developed from the following considerations. According to 
equations (73) and (9*0, 
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for V = — . 

2 

that Q = V 


when V = — 
2 

computation 


Therefore, 


and 


Dll '=-l-2oS 


■ D 12 ' ~ D 22 ' = 1 



(06) 


From these results and equations (92), it is easily seen 
consequently. 


k = | (1 - 3c8) 


(07) 


. For this value of Poisson's ratio, a straightforward 
will yield 


c = 


a 


1-1 
31 

_ - *> + 1 
1 - 1 


5 + = - .)TL- 1 
° sfr+l 

8 = ^ ^ 

" (l + \fX) 3 




(08) 


cB = (V * : (09) 

3(l +n/i) 2 


k = | (1 - 3c5) 



(cio) 
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On the other hand, according to equations (9*0* 


9g 0 a 2 

A% 0 


(cn) 


■where Is the critical compressive stress. Consequently, 

Ci 


, 1 

3c 2 (l + ^) £ 

and the critical buckling load P is given by 


( 012 ) 


P = 


°cr* h 


A. 1 4 

c 2 (l+\fx) 2 


from equation (99) * This is the desired result from which Von Kerman's 
eq.ua tion follows. 

It has been pointed out in the third section of the ANALYSIS that 
the neutral line of the beam and the neutral surface of the plate will 

coincide if and only If V = ^. This can be seen quite simply from 

the following discussion. If V = ^, then z Q = according to 

equations (^9) • It follows from equations (c8) that 


t + = ffo \[X - 1 

° h \/I + l 


(C13) 


The distance of the neutral surface from the lower surface is given by 


h + z = 
2 + Z ° 


h 


= h 




\pK + 1 + \[e 


(C14) 
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This is the same result as that obtained in the direct analysis of the 
buckling of a beam.. (See reference 1, p. 158.) 

Conversely, assume that the two positions of the neutral line 
agree; that is. 


t + - - VT - 1 

0 \[x + 1 


(C15) 


where £ Q + is computed from equation (Bl6) . Since equation (C15) 

must be an identity in X, it must be valid for any particular value 
of X. For computational convenience, consider the case when X = t. 
Then 


\f\. - 1 _ 1 

Vx + 1 3 


c 


2 

19 - 12V - 4v 2 




(C16) 


a = 3(5 Hv) 


where u 


8v 2 + 12V 


23 • Consequently, 


+ 3 - 9(5 - 4V)2 

_ 1 

3 


(C17) 


and subsequent simplification leads to 


u = -5(5 " ^V) 


(Cl8) 



NACA TOT Wo. 1530 


79 


Evaluation of u in terms of V yields 


(2V - l ) 2 = 0 


(C19) 


or 


V = 


1 

2 


(C20) 


Therefore the two methods of determining the neutral line will agree 
if and only if v = and the transition from plate to team will 

"be valid only under this condition. 


Buckling of a Cruciform Section 

The transcendental equation ( 138 ) must he solved to determine 
the critical stress. A straightforward computation will show that 


S 2 


m 2 v% 2 

B22 ' 




2 

n 


mW 

B22' 


^°22 ' 


m 


AD, 


22 



■s 




(C21) 


Since 


i 2 = T,e 


2m 2 iPsPv. 


12 


22 


(C22) 


Equation ( 138 ) can he written entirely, in terms of tj and the 
parameters E, Q, > and D 22 ' ’ ^ or values of m and E, 

the remaining parameters in the transcendental equation are functions 
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of X alone. Hub tj can be determined as a function of X by 
solving this equation by the usual iterative procedures. The 
quantity k can then -be found as a function of X by solving the 
second of equations (C2l) . These results are independent of the 
stress -strain lav of the plate material. On the other hand, 

k - (l - v 2 ) (023) 

« 2 h 2 E Q 

that is, for a given material k can be written as a function of X 
once ct q is known as a function of X ■ 

It is a little more convenient to write the solution in terms of 
the quantity k rather than _k. This parameter was defined in 
equation (117) • Curves for k as a function of X for given values 
of m and E can be obtained in the following way. The quantity k 
can be computed, as previously described, by solving equations ( 138 ) 
and (C21); E can then be found from equation (117) • -These results 
are independent of the stress-strain law and_cansequently hold for 
any rectangular plate. On the other hand, k can also be computed 
as a function of X for a given stress-strain lav according to 
equation ( 117 ) • The resulting function will depend on the 
parameter b 2 /h 2 . Once the side ratio E, the width-to- thickness 
ratio b/h, and the wave form m have been selected, the value of k 
corresponding to the buckling stress can be obtained by finding the 
Intersection of the two corresponding curves of E against X 
computed as Just outlined. 

It has been pointed out in the third section of the ANALYSIS 
that the lowest value of k, .and hence the lowest critical stress, 
will be attained for m = 1. This can be seen easily from the 
following considerations . Figure 9 showed that for m = 1 and a 
fixed value of X, k increases for increasing R. With a simple 
change of variables, the solution for any value of m can be obtained 
from the solutions for m = 1. Let 



E ' = mR 

E 



(C24) 
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Then equations (138) and (C 21 ), written in terms of k* and E' will 
"be precisely the same as those for k and B when m = 1 . Thus the 
curves of figure 9 can he used for any m. provided E is replaced 

by E' and k by k f . For m>l then, k > k 1 and E < E * • In 
other words when m. ^ 1 , the value of k for a given value of E 
and X will be larger than the corresponding value of k for m = 1. 
Since the dashed curves are monotonic, increasing functions of X, 
this implies that the critical value of k will be lowest for m = 1* 
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TABLE I 


YALOES OF D-^’, AND D £2 ’ 

FOR V = 0.32 


A 

D n* 

d 12 ’ 

»22' 

1.0 

1.0000 

1.0000 

1.0000 

1.2 

•93030 

1.0119 

.99680 

11 

.87II7 

1.0269 

.99421 

1.6 

.828I7 

1.0368 

.99212 

1.8 

.78973 

1.0151 

.99031 

2.0 

.75655 

1.0522 

.98882 

2-5 

.69086 

1.0662 

.98580 

3-0 

.6H78 

1.0768 

.98355 

3-5 

.603I5 

I.0850 

.98179 

1.0 

•57255 

1.0916 

.98037 

1.5 

.5I698 

1.0971 

.97920 

5-0 

.525H 

1.1017 

.97821 

6 .0 

.I910I 

1.1091 

.97663 

7-0 

.I6I59 

1.1117 

.97511 

8.0 

.II365 

1.1192 

.97^15 

9-0 

.I2652 

1.1229 

•97367 

10.0 

.I1226 

1-1259 

.97301 
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Figure 4. - The quantity k/ (l - v 2 ) 
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as a 
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ction of X for two beams having the same shape 
; values of v . 
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Figure 5.- Stress -strain curve for e = -^-<x + Acr^ , E Q = 10,667,000, 

E o 

log A = -42.873, and 0 = 8.6127. 
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(a) m = 1. 

Figure 6. - Curves of k as a function of A for different values of m. 
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Figure 6.- Continued. 









NACA TN No. 1530 



(c) m = 3. 


Figure 6.- Continued. 
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(d ) m = 4. 


Figure 6.- Concluded. 
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Figure 7.- Comparison of values of k for values of 1/R with 
v = 0.32 as determined by different theories. 
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Figure 8.- Comparison of values of k for values of 1/R with 

v = 1/2 as determined by the present method and that of Ilyushin. 
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Figure 9. - Results of solving for k as a function of 

equations (117) and (138). 


A by using 
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Figure 11.- Curves of k obtained by the energy method and exact 
solution for v = 0.32 and various values of R. 
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Figure 12.- Results of applying energy method to cases of figure 10. 






